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Multi Channel Grating Design 



Field of the invention 

The present invention relates broadly to a multi-channel grating design method land to 
multi-channel grating structures fabricated utilising the multi-channel gratmg design method. 

Background of the invention 

Multi-channel grating structures are typically written into photosensitive waveguides. 
The grating structure comprises a refractive index profile induced in the photosensitive 
waveguide, which in turn determines the optical characteristics such as the reflection, 
transmission, and group delay characteristics of the resulting gratmg structure. 

The refractive index profile for a multi-chaimel grating structure is induced by applying 
a sampling function to a given single-channel grating profile, the single-channel grating profile 
typically being a periodic profile achieved by exposing the photosensitive waveguide to an 
appropriate light beam emerging from e.g. a suitable phasemask. 

In a vast majority of previously reported work on multi-channel gratings, a so-called 
Sinc-sampled design has be^ used. For the Sinc-sampling approach, an i\r-channel gratmg 
design can be obtained by a direct in-phase summation of ^ idmtical seeding gratings [with 
k(z) - amplitude grating amplitude, 0 (z) - grating phase] equally spaced in the frequency space: 



Gsinc =Xr-i^^^t(2/ -iV- l)Afe/2] = " 2ra)/2],sinc(x) s sin(;c)/x , 

and Ak is the channel spacing. 

This design will be referred to as "in-phase" grating design herein after by the applicant. 

An example of this design is shown in Figures 9 (c) and (d) with corresponding spectral 
characteristics as shown in Figures 9 (a) and (b). The maxhnum value of the refiractive index 
change required to implCTient this multi-channel grailmg design is given bya siiiiplei^xpression: 




(1) 



where 



(2) 
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where An, is the maximum refractive index change required for the single seeding grating. 
Since any photosensitive fiber iised to fabricate Bragg gratings has material limits of the 
maximum achievable photoinduced refractive index change An^ this represents a limitation on 
the maximum number of channels that can be recorded in a given fiber. Thus it is highly 
desirable to reduce a required An^^ as much as possible. 

At least prefCTred embodiments of the present invration seek to provide an alternative 
multi-channel grating design in which the maximum refractive index change required as a 
function of the number of channels is reduced without deterioration of the grating spectral 
characteristics when compared with the prior art grating designs discussed above. 

Summary of the invention 

in accordance with a first aspect of the present invention there is provided a method of 
calculating a sampling function for fabricating a ^-channel grating, the method comprising the 
steps of forming a summation of iV^ periodic seeding functions each describing a refractive index 
variation, wherein each periodic function includes a phase shift value with respect to the other 
functions, and wherein at least one of the phase shift values is non-z^o. 

Accordingly, the present iuvention can provide a "de-phased" grating design rather than 
the "in-phase** prior art grating design. 

Preferably, the summation of tihe N periodic functions comprises a Fourier analysis. The 
result of the Fourier analysis may be expressed as: 

where Q is the ampUtude and \\f is the phase of the samphng function. 

The method may further comprise the step of determining a set of the phase shift values 
for which a maximum of the sampling function ampUtude Q is minimised. 

Alternatively, the method further may comprise the step of determining a set of the 
phase shift values for which a maximum difference between a minimum and a maximum of the 
samplmg function amplitude Q is minimised. 
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Alternatively, the method may further comprise the step of determining a set of the 
phase shift values for which a mean-square-deviation in the sampUng function amplitude Q is 
minimised. 

The step of determining the set of phase shift values may comprise direct scanning 
through all combinations or conducting a variational analysis, or using other forms of extremum 
search numerical techniques, or a simulated aimealing Monte Carlo approach. 

The grating may be multi-dimensional, wherein the periodic seeding functions are multi- 
dimensional. 

In accordance with a second aspect of the present invention, there is provided a method 
for fabricating a multi-channel grating comprising the step of calculating a sampling function in 
accordance with a method as defined in the first aspect. 

The multi-channel grating may e.g. be fiibricated utilising photo-induced refiractive index 
changes in a photosensitive waveguide material, etching techniques, or epitaxial techniques, or 
a developing technique such as a photo polymerisation process. 

In accordance with a third aspect of the present invention, there is provided a multi- 
channel grating structure fabricated utiUsing a method of fabrication as defined in the second 
aspect. 

Brief description of the drawings 

Preferred forms of the present invention will now be described with reference to the 
accompanying drawings. 

Figure 1 shows a four-channel grating design using the sampling function with 
minimised maximum amplitude embodying the present invention. 

Figure 2 shows the maximimi value (maximum minimisation) of the sampling function 
ampUtude for different nimiber of channels, (a) Comparison with the worst case (in-phase 
sampling function design), (b) Comparison with the asymptotic theoretical value ^In . 

Figure 3 shows a 16-channel grating design using maximum minimisation ^proach 
mbodying the present invention. 
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Figure 4 shows the maximum value (difference minimisation) of the sampling function 
amplitude for different number of chamiels. 

Figure S shows a 16-channel grating design using the difierence minimization s^proach 
embodying the present invention. 

Figure 6 shows a 17-channel grating design using the variational minimization approach 
embodying the present invention. 

Figure 7 shows the mflyiyni'Tn value (variational approach) of the sampling function 
amplitude for a difiG^ent number of channels. 

Figure 8 shows an experimental set up for writing a multi-chamel grating structure of a 
multi-channel grating design mibodying the present invention. 

Figure 9 shows a prior art four-channel grating design using an in-phase-samplmg 
function. 



Detailed description of the ennbodiments 

The preferred embodiment described provides a multi-channel grating design in which 
the maximum refiactive index change is less than directly proportional to tiie number of 
channels N, and which avoids the presence of un-wanted side bands in the spectral 
characteristics of the resulting grating, thereby improving on prior art multi-channel grating 
designs. 

In the preferred embodiment, a sampling fimction which periodically modulates the 
amplitude of a given single-channel grating (seeding grating) is utilised, similar to prior art 
multi-chaimel grating designs. However, in addition to the periodic modulation of the 
amplitude of the seeding grating, a periodic modulation of the phase of the seeding grating is 
also introduced. Accordingly, the resulting design function in the preferred embodiment may be 
expressed as: 

where the phase of the sampling function t|/=\|;(2) and the sampling amplitude Q=Q(z) are 
given by: 
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g2 ^ 4 2]cos(a, - jcos(/i; Afe / 2 + ) cosCn^ Afe / 2 + ), and 

/.pa] 



^(z) = tan 



-1 



J cosa^cos(n^Afe/2 + /?^j 



is even 



or 



(AM)/2 

e^(2) = 4 2]cosa,cos(n;Afa/2 + P;) + 

4 2]cos(a; -a^)cos(n;Afe/2 + p;)cos(n^Afe/2 + p^)+l, 



and 



^/(z) = tan"* 



^l^^^^'^sina, cos(/i;Afe/2 + P;) 



JVisodd, 



^|^j"^^'^cosa| cos(/i,M2r/2 + p, + 1) 
where s2/~^-l and 71^ s2/7"-iV-l. 

In the above expressions for Q(z) and y/(z) we use notations 
a s (^^ + ^i*r+i-/)/2, s (^^ -^i^^,_^)/2 and set $>(jv+i)/2 = 0 for odd number of channels. Now for 
any given there will be a set of {a^},{P;} (orequivalently aset of (cj),} ) which optimizes ttie 
sampling function amplitude g(z) along the grating structure according to a selected 
optimization approach. By direct calculations it is straightforward to show that 

Anitk 



(4) 



for any choice of ai and pi. This expression, in turn, leads to an asymptotic formula for the 
minimum possible An^ corresponding to an "ideal" situation when Qiz)-4N and only the 
grating phase is nontrivially modulated (by addition of an approprietly chosen phase of the 
sampling function v^(z)): 



(5) 



We note, that, in practice, the limit Q(z) =->/JV can be reached only approximately. For 
example, for the maximum minimization approach, one should solve a minimax problem and 
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find Q^(z;aj""\fil'^'^)foT which max,{8^(2;ar>.pr) } = mm(.,.p,,max,{e(z;a„p,)}. 
To find the optimal set for a relatively small number of channels one may use direct 

numerical scanning through all possible combinations of the dephasing angles. An example of 
the optimized "out-of-phase" design for a four-channel grating is shown in Fig. l(a)-(d). For 
this design a, = 0.5759;r,a2 = = = 0 . This design reduces the maximum An^^ of the in- 
phase design by aknost 40%. 

For A^»l location of the minimizing set (a/**'"\>9/''^'^)is a nontrivial exercise. Even 

rough direct scanning through all possible sets of angles (followed by efficient niunerical 
minimum search routines) quickly becomes numerically inefficient. To solve the optimization 
problem for large N we use, in an example OTibodiment the so-called simulated annealing 
method - a Monte Carlo £^proach for minimization of multi-variable fimctions. This statistical 
method samples the search space in such a way that there is a high probability of finding an 
optimal or a near-optimal solution in a reasonable time. The term ''simulated annealing'" is 
derived fix>m the analogy to physical process of heating and then slowly cooling a substance to 
obtain a crystalline structure. 

To start, the system state is initialized. A new configuration is constructed by imposing a 
random displacement. If the energy of the new state is lower than that of the previous one, the 
change is accepted and the system is updated. If the energy is greater, the new configuration is 
accepted with some probability. This procedure allows the system to move consistently towards 
lower energy states, yet still jump out of local minimal due to tibe probabilistic accq>tance of 
some upward moves. Results of the application of simulated annealing to the maximum 
minimization approach are sunmiarised in Figure 2. 

Figure 2 shows that Eq. (5) gives a reasonable estimate for An^ (20% accuracy or 

better) for N > 6. The growth factor to ^^N represents a significant improvement in 
comparison to linear dependence (almost a factor of 4 reduction for 16-channel gratings, 8 
for 64-channel ones and so on). An example of the optimal 16-channel grating design is shown 
inFig.3(a)-(d). 

Figures 1 and 3 demonstrate implementation of Q^(z) reduction strategy without 
trying to avoid touching the zero level at some z. Zeros in the fibre Bragg grating (FBG) 
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amplitude may lead to the increased phase errors (appearance of phase jumps) and should be 
avoided. Thus, arguably, the best minimization strategy could be not in reducing the maximum 
value of Qiz) along z, but in minimizing the maximum difference betweai the maximum and 
TniniTnum of the sampling function amplitude. Mathematically this may be formulated as 



finding Q^(z;aj^''\j3}^''') for which m^,^g^iz;a^^p^''')^^ 
min^„,,p^)[max,{e(2;a;,p;)}-min,{e(z;a;,p,)}]. 



This ^proach may be implemented by using the same simulated annealing algorithm, 
and the corresponding results are given in Fig. 4. An example of this version of optimal 
(difference minimisation) "out-of-phase'* design for a 16-channel grating is shown in Fig. 5(a)- 
(d). The total minimization of Ahj^ is slightly worse than in Fig. 3, but the zeros in the grating 
amplitude are avoided. 

It is noted that for both the- maximum niinimisation and the difference minimisation 
^proaches, points with the best optunisation quality may be used to get reasonable (but not 
exactly the best) optimisation for some higho: channel numbers. For example, the point N=0 
vnfh the particularly good optimisation quality provides us with an effortless optimisation for 
i^=81 number of channels. 

Another embodiment will now be described, in which optimization by the functional 
minimisation (variational approach) is utiUsed. The key property of this embodiment is that it 
relies on estimate of some integral fimctional rather than time-consuming numerical scanning in 
z. Quantitatively, proximity of g(z) to the theoretical limit ^In can be characterised by mean- 
square-deviation, 




(6) 



where {f(z)) ^f-tskl2n^ f{z)dz Ideal optimisation of Q{z) corresponds to the 

achievement of the average Q = ViV and the zero meau-square-deviation fi:om this average 
value. Using expression (4) and assuming Q « 4n , one can show that 



(7) 



where 
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1 N N N4l2P\ 

/si Bsl mat 



cos(a„ - a 



.^i-p\)^^i^i +P/ "P^ +|J„ -P„,|,_^,) . 



For finding minima of Aj2 the most efficient strategy is again the use of the simulated 
annealing method. 

The major advantage of the optimisation based on the functional minimisation compared 
with direct scanning is the speed: integration over z is carried out analytically which saves lots 
of computer time. However, the quality of variational optimisation itself is usually not as good. 
That said, for odd nimiber of chaonels variational optimisation leads to a sampling function 
without zeros in amplitude (similar to diflfeence minimisation sqpproach). An example for a 17- 
diannel grating design (variational ^proach) is shown in Fig. 6(a)-(d). For even number of 
channels variational optimisation leads to a sampling function with zeros in ampUtude (similar 
to maximum minimisation approach)..Fig. 7 shows-the maximum value (variational approach) 
of the sampling function amplitude for a different number of channels. 

The advantages of the dephasing optimisation ^proaches of the embodiments described 
include that they give a very significant reduction for the maximum required An^ . Moreover, 
they are conceptually simple and relatively easy to obtain. Furthermore, the limitation of the 
number of dephased periodic functions to N for a iV-channel grating design avoids the presence 
of im-desired side bands, which can deteriorate the quality of spectral characteristics of the 
grating. The q}plicants have recognised that the limiting of the number of dephased periodic 
functions to while allowing variations in the amplitude of the resulting sampling function can 
enable design of ^-channel gratings of improved quality when compared with prior art designs. 

The implementation of the multi-chaimel grating design of the preferred embodiment in 
a grating structure requires grating writing apparatus with high spatial resolution to be utilised. 
Therefore, in a grating writing apparatus relying on photoinduced refractive index changes, the 
apparatus preferably comprises a beam focusing means to reduce the size of the beam in the 
core of the photosensitive waveguide. 

Figure 7 shows an example experimental set up SO for writing a multiTchannel grating 
52 into an optical fibre 54. The experimental set up 50 comprises an inteiferometer 56 which 
includes a first acousto-optic modulation 58 being operated under an acousto-optic wave at a 
first fiequency Qi, as indicated by arrow 14. An incoming li^t beam 60 is incident on the first 
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acousto-optic modulator 58 at a first order Bragg angle. The operating conditions of the 
acousto-optic modulator 58 are chosen such that the modulator 58 is under driven, whereby 
approximately 50% of the incoming beam 60 is diffracted into a first order beam 62, and 50% 
passing through the acousto-optic modulator 58 as an un-dif£racted beam 64. The un-difi&acted 
beam 64 is incident on a second acousto-optic modulator 66 of the interferometer 56 at a first 
ordCT Bragg angle, whoreas the beam 62 is not. Accordingly, the beam 62 passes tfarougih the 
second acousto-optic modulator 66 without any significant loss. 

The second acousto-optic modulator 66 is operated under an acousto-optic wave at a 
frequency Q2» which propagates in a direction opposing the direction of the acousto-optic wave 
in the first modulator 58 as indicated by arrow 68. After the second acousto-optic modulator 66 
the first order dififracted beam 70 and the beam 62 are firequency shifted in the same direction 
(e.g. higher fi-equency), but by different amounts i.e. £l\ v Q2- 

The beams 62, 70 are then brought to interfere utilising an optical lens 72, and the 
resulting interference pattern (at numeral 74) induces refiiactive index changes in the 
photosensitive optical fibre 54, whereby a refiactive index profile, i.e. grating structure 52, is 
induced in the optical fibre 54. 

In Figure 7, the optical fibre 54 is translated along the interferometer at a speed v, as 
indicated by arrow 74. 

It will be appreciated by a person skilled in the art that the experimental set up 50 shown 
in Figure 5 can be utilised to write a multi-channel grating structure of a multi-channel grating 
design embodying the present invention through suitable control of the first and second acousto- 
optic modulators 58, 66, in conjunction with a suitable control of the speed v at which the 
optical fibre 54 is translated along the interferometer 56 at any particular time. The high spatial 
resolution required to implement the multi-channel design of the preferred embodiment is 
achieved in the set up shown in Figure 6 by utilising the optical lens 72, with the practical limit 
of the beam size in the focal plane preferably being of the order of the waveguide core size. 

It will be appreciated by a person skilled in the art that numerous variations and/or 
modifications may be made to the present invention as shown in the specific embodiments 
without departing fi-om the spirit of scope of the invCTition as broadly described. The present 
embodiments are, therefore, to be considered in all respects to be illustrative and not restrictive. 
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For example, multi-channel gratings can be fabricated on ihe basis of the multi-channel 
grating design of the present invention using various known grating writing techniques, 
including one or more of the group of photo-induced refractive index variation in photo 
sensitive waveguide materials, etching techniques including etching techniques utilising a 
phasemask, and epitaxial techniques. Furtiiermore, while the preferred embodunents have been 
described in the context of l-dimrasional Bragg gratings, the present invention does extend to 
multi-dimensional multi-channel gratings. Such gratings have applications e.g. as photonic 
bandgap structures. 
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1 . A method of calculating a sampling function for fabricating a AT-channel grating, 
the method comprising the steps of: 

- forming a summation of M periodic seeding functions each describing a refractive 
index variation, wherein each periodic function includes a phase shift value with respect to the 
other functions, and wherein at least one phase shift value is non-zero. 

2. A method as claimed in claim 1, wherein the summation of die N periodic 
functions comprises a Fourier analysis. 

3. A method as claimed in claim 2, wherein the result of the Fourier analysis is 
expressed as: 



4. A method as claimed in any one of claims 1 to 3, wherein the method furth^ 
comprises the step of detamining a set of the phase shift values for which a maximum value of 
the sampling function amplitude is minirnised. 

5. A method as claimed in any one of claims 1 to 3, wherein the method further 
comprises the step of determining a set of the phase shift values for which a maximum 
difference between a maximum and minimum value of the sampling function amplitude is 
minimised. 

6. A method as claimed in any one of claims 1 to 3, wherein the method further 
comprises the step of detennining a set of the phase shift values for which a meau-square- 
deviation in the sampling function is minimised, 

7. A method as claimed in any one of claims 4 to 6, wherein the step of determining 
the set of phase shift values comprises direct scannmg through all combinations, or conducting 
a variational analysis, or using other forms of extremum search numerical techniques, or a 
simulated annealing Monte Carlo approach. 

8. A method as claimed in any one of the preceding claims, wherein the •grating is 
multi-dim^isional, and wherein the p^odic seeding functions are multi*dimensional. 



1=1 
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9. A method for fabricating a multi-channel grating comprising the step of 
calculating a sampling function in accordance with a method as claimed in any one of claimis 1 
to 8. 

10. A method as claimed in claim 9, wherein the multi-channel grating is fabricated 
utilising photo-induced refractive index changes in a photosensitive waveguide material. 

11. A method as claimed in claim 9» wherein the multi-channel grating is fabricated 
utilising etching techniques. 

12. A metibiod as claimed in claim 9» wherein the multi-channel grating is fabricated 
utilising epitaxial techniques. 

13. A method as claimed in claim 9, wherein the multi-channel grating is fabricated 
utilising a developing technique. 

14. A method as claimed in claim 13, wherein the developing technique comprises a 
photo polymerisation process. 

15. A multi-channel grating stracture &bricated utilising a method of fabrication as 
claimed in any one of claims 9 to 14. 
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